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ABSTRACT 

Using a self-consistent atmosphere code, we construct a new model of the atmosphere of the transiting 
extrasolar giant planet HD 209458b to investigate the disparity between the observed strength of the 
sodium absorption feature at 589 nm and the predictions of previous models. For the atmospheric 
temperature-pressure profile we derive, silicate and iron clouds reside at a pressure of several mbar in 
the planet's atmosphere. These clouds have significant vertical extent and optical depth due to our slant 
viewing geometry and lead to increased absorption in bands directly adjacent to the sodium line core. 
Using a non-LTE sodium ionization model that includes photoionization by stellar UV flux, coUisional 
processes with H2, and radiative recombination, we show that the ionization depth in the planet's 
atmosphere reaches ~l/2 mbar at the day/night terminator. Ionization leads to a slight weakening of 
the sodium feature. We present our baseline model, including ionization and clouds, which falls near 
the observational error bars. The sensitivity of our conclusions to the derived atmospheric temperature- 
pressure profile is discussed. 

Subject headings: planetary systems — radiative transfer — stars:individual(HD 209458) 



1. INTRODUCTION 

The planet in the system HD 209458 has been an object 
of intense study since the first observations of its transit 
(Henry et al. 2000; Charbonneau et al. 2000). These ob- 
servations led to a determination of the planet's radius and 
physical mass, confirming that it must be a gas giant. At 
wavelengths where the atmospheric opacity is high, the 
planet blocks more light, yielding a deeper transit light 
curve that makes the planet appear physically larger. The 
location and strength of these absorption features serve 
as diagnostics of the temperature and chemistry of the 
planet's atmosphere. Seager & Sasselov (2000) predicted 
variations in the transmission spectrum of HD 209458b 
due to neutral Na and K and to singly-ionized He. Brown 
(2001) explored effects such as ionization and winds and 
predicted large variations due to H2O, CO, and CH4 in 
the infrared. Hubbard et al. (2001) explored a variety of 
physical effects that could be found in transits, including 
refraction and the angular redistribution of photons due 
to Rayleigh scattering. These were found to be minimal 
for HD 209458b; Hubbard et al. (2001) went on to derive 
the planet radius as a function of wavelength from 300 to 
2500 nm. 

Confirmation that absorption features due to the opac- 
ity of gaseous species could be observed with current tech- 
nology was obtained by Brown et al. (2001). They ob- 
served the transit of HD 209458b from 581 to 638 nm with 
the STIS instrument on HST and obtained photometric 
accuracy near 100 micro-magnitudes. This spectral region 
encompasses the Na-D doublet at 589 nm. Using these 
data, Charbonneau et al. (2002) found that the transit was 



deeper by (2.32 ± 0.57) x 10~^ in a narrow band centered 
on the star's Na-D lines than in adjoining bands at shorter 
and longer wavelengths. Importantly, the observed differ- 
ence in transit depth in and out of the Na-D line is smaller 
than previously predicted (Charbonneau et al. 2002). The 
nominal model of Hubbard et al. (2001) predicted a differ- 
ence between these two bands of 4.7 x 10"''. 

Charbonneau et al. (2002), using the models of Brown 

(2001) , performed a parameter study to investigate a vari- 
ety of possible reasons for this difference. These included 
a bulk underabundance of Na, the sequestering of atomic 
Na in condensates and/or molecules, clouds very high in 
the planet's atmosphere, and photoionization of atomic Na 
due to the UV flux from the parent star. Barman et al. 

(2002) suggested that Na is not in local thermodynamic 
equilibrium (LTE). They speculated this leads to a core 
reversal of the 589-nm absorption doublet in the planet's 
emission spectrum, but they did not calculate its effect on 
a transit light curve. 

Here we explore the effects of cloud formation and ion- 
ization of Na in HD 209458b's atmosphere. We flnd that 
cloud opacity is the dominant effect that determines the 
depth of the Na-D feature, but ionization of Na also leads 
to a non- negligible weakening of the Na-D doublet. 

In our analysis, following Charbonneau et al. (2002), 
we define a narrow "in feature" wavelength band as the 
wavelength range 588.7-589.9 nm, and an "out of feature" 
band as the combination of the 518.8-588.7 and 589.9-596.8 
nm wavelength ranges. The "in feature" band includes 
the wavelength extent of the easy-to-observe HD 209458A 
(stellar) Na absorption doublet. When quoting a depth 
for the Na-D feature, we simply mean the difference in 
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transit depth (the fractional change in flux) between the 
"in feature" and "out of feature" bands at midtransit. The 
fractional change in flux at midtransit across the entire 581 
to 638 nm wavelength range is -0.0164, meaning 1.64% of 
the star's light is blocked. 

2. TEMPERATURE-PRESSURE PROFILE 

The temperature-pressure ( T-P) profile for a close-in ex- 
trasolar giant planet is flatter (more isothermal) than that 
of an isolated object with the same effective temperature 
because of heating of the outer atmosphere by the parent 
star. In the present work, we have employed TLUSTY 
(Hubeny & Lanz 1995), a self-consistent atmosphere code 
that uses the Discontinuous Finite Element (DFE) method 
for radiative transfer (Hubeny & Lanz 1995; Castor et al. 
1992). The atmosphere model is plane parallel. Radia- 
tive equilibrium is achieved by adjusting the atmospheric 
temperature-density structure iteratively via the Accel- 
erated Lambda Iteration (ALI) method. The important 
physical quantities, including net flux, temperature, and 
density in every atmospheric depth zone, are converged 
to better than 0.1%, and total energy is conserved to this 
same precision. For the deep convective layers a mixing- 
length prescription is used with a mixing-length parameter 
of one scale height. 

For a close-in planet such as HD 209458b, an inner 
boundary flux corresponding to that of an isolated object 
with an effective temperature of ~500 K or less results 
in an outer atmospheric structure that is independent of 
this inner boundary. The outer structure is determined 
by the stellar flux, the surface gravity, the atmospheric 
compositions, and the opacities (see Sudarsky et al. 2002; 
Burrows et al. 2001). The metallicity of the parent star 
has been shown to be consistent with solar abundances 
(Mazeh et al. 2000), and we use these abundances when 
modeling the planet's atmosphere. 

The T-P profile of HD 209458b's atmosphere that we 
derive is shown in Fig. 1. We use a surface gravity of 
980 cm/s^ for the planet. We find that HD 209458b is 
a "Class V" planet as defined in Sudarsky et al. (2000) 
and Sudarsky et al. (2002), meaning that silicate and iron 
clouds occur at altitudes in the planet's atmosphere rele- 
vant to spectroscopic observations. We model the location, 
distribution, and size of cloud particles using the model of 
Cooper et al. (2003). As seen in Fig. 1, the cloud bases of 
forsterite (Mg2Si04) and iron (Fe) are at the several mbar 
level. The effect of these clouds on the transit is discussed 
in Section 4. 

In calculating the T-P profile, we assume that incoming 
stellar energy is absorbed and redistributed evenly on the 
day side, but is not redistributed onto the night side. In a 
recent paper. Showman & Guillot (2002) find that redis- 
tribution of energy around the entire planet may be ineffi- 
cient, leading to a ~500 K effective temperature difference 
between the planet's day and night sides. The transit ob- 
servation is sensitive to the planet's atmosphere on the 
limb, which could have a T-P profile slightly cooler than 
we predict here. (Note that the atmospheric circulation 
code of Showman & Guillot (2002) is not able to calcu- 
late T-P profiles for the planet.) Based on this possibly 
large temperature difference, Guillot & Showman (2002) 
make the point that if HD 209458b does have a 500 K 



effective temperature contrast between the night and day 
sides, condensates could form on the night side and tie 
up Na, thereby weakening the Na-D absorption feature. 
Note, however, that the "cold" T-P profile used in Hub- 
bard et al. (2001), with a temperature 400 K cooler at 1 
mbar than seen in Fig. 1, still has a 589-nm Na feature 
depth of 3.2 x 10""*, deeper than observed. 

Our T-P profile differs significantly from that calcu- 
lated in the recent paper on HD 209458b by Barman et al. 
(2002). We find a temperature of 1280 K at 1 mbar and 
950 K at 10 /Ltbar. Examination of Fig. 1 in Barman et al. 
(2002) shows temperatures of ~1790 K at 1 mbar and 1630 
K at 10 /ibar. The reasons for these large differences are 
unclear, but could be due to differing opacities, stellar flux, 
or planet surface gravity. Our T-P profile is closer to those 
for HD 209458b discussed by Seager & Sasselov (2000) and 
Richardson et al. (2002). These differences are significant, 
and we are led to very different conclusions concerning 
the physical processes of importance for the Na-D transit 
depth. As with late L and T dwarfs, our HD 209459b at- 
mosphere is free of metals such as Mg, Al, Ca, Fe, and Si, 
which have been incorporated into condensates. Barman 
et al. (2002) find an atmosphere that is not cool enough to 
allow for the formation of these condensates, which leaves 
their atmosphere rich in heavy metals in atomic form. Due 
to the high UV opacity of free metals that remain in the 
atmosphere and absorb stellar radiation. Barman et al. 
(2002) conclude that ionization of Na at pressm-cs greatc;r 
than 1 /ibar is negligible. Our model temperatures allow 
deeper ionization of Na in the atmosphere of the planet. 

3. TRANSIT CALCULATIONS 

The methods we use for calculating the effects of ab- 
sorption by gaseous species, absorption and scattering by 
condensates, Raylcigh scattering, and refraction are essen- 
tially identical to those described in Hubbard et al. (2001), 
although we remove the simplification that the opacity ver- 
sus pressure profile is the same around the entire planet. 

We calculate the slant optical depth through the planet's 
atmosphere, r, from the expression 
^_ r r"dr"a[r")N[r") j"^ r"dr"a{r")N(r") 

^ ~ J-oo (r"2 - r2)^/^ Jr (r"2 - r2)^/^ ' 

where r is the impact parameter of a light ray, r" is the 
magnitude of the three-dimensional vector position of an 
atmospheric point from the planet's center, a is the ab- 
sorption cross-section per molecule, and N is the number 
density of molecules in the atmosphere. Two terms in the 
T expression are needed because cr(r") is different on the 
day and night sides. 

As before, we use the planetary and stellar properties 
from Brown et al. (2001). Fig. 2 shows, for a cloudless 
atmosphere with no ionization, the radius as a function 
of wavelength over the range covered by the Brown et al. 
(2001) and Charbonneau et al. (2002) observations. The 
auxiliary axis of Fig. 2 shows the pressure at the termi- 
nator to which this radius corresponds; the figure thus 
indicates the pressure range probed by the transit of HD 
209458b. 

4. CLOUD MODEL 

4.1. Location, Distribution, and Opacity 
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The T-P profile we derive is warmer than the a,d hoc 
profiles studied in Hubbard et al. (2001). Because of this, 
we find forsterite and iron clouds higher in the planet's 
atmosphere. If we examine the forsterite cloud in Fig. 1, 
for the T-P profile without cloud opacity, the cloud base 
is at 5.4 mbar. The opacity due to clouds slightly warms 
the planet's atmosphere, especially at the location of the 
clouds, moving the cloud base farther out, to 1.8 mbar. 
Predicting the exact location of the cloud base is a non- 
trivial problem, and in a fully self- consistent calculation 
of cloud location, one would need to iterate many times 
successively between the cloud code and atmosphere code. 
Since cloud opacity is a strong function of particle size, and 
the predicted particle sizes are altered due to a change in 
the cloud base and as function of pressure within the cloud 
itself (as described below), this step is challenging. Since 
uncertainties remain in the modeling of clouds, and the 
effects of redistribution of energy on the T-P, we do not 
attempt this fully self-consistent treatment. Instead, our 
T-P profile shown in Fig. 1 incorporates the opacity of the 
cloud with the location of the base given by the point of 
intersection with the condensation curve, and with a cloud 
particle size given by the modal size for the entire cloud. 
The particle number density falls of exponentially for one 
scale height. 

The cloud model of Cooper et al. (2003) predicts par- 
ticle radii in this radiative region for the iron cloud that 
range from tens of /im at the base to less than a /im at 
the cloud tops. In the forsterite cloud, all particles are 
0.05 to 0.5 nm in size. The cross-sections are calculated 
as in Sudarsky et al. (2002). For these submicron parti- 
cles, scattering dominates over absorption, and the opac- 
ity of the forsterite cloud is significantly larger than that 
of the iron cloud. The Cooper et al. (2003) model pre- 
dicts a distribution of cloud particle sizes and densities as 
a function of pressure in the atmosphere. Therefore, we 
calculate the cross-sections for the cloud as a function of 
pressure at many pressures within the cloud and put no 
artificial cap on the cloud top. While understanding the 
vertical distribution of cloud material is a well-known area 
of importance for brown dwarf atmospheres, it is of even 
greater importance for the transit spectroscopy of EGP at- 
mospheres, because cloud opacity can be so great for our 
slant viewing geometry. 

For HD 209458b, we find a forsterite cloud with sig- 
nificant vertical extent, with a large slant optical depth (r 
~500 at the base) that does not fall to r ~ 2/3 until 1.9 gas 
scale heights above the cloud base. Even though 99% of 
the cloud's mass is found within 1 gas scale height, a non- 
negligible population of small, highly scattering cloud par- 
ticles continues above this level, giving rise to a forsterite 
cloud with a significant slant optical depth up to ~1.9 
scale heights. Since our predicted forsterite cloud is in a 
region of interest for HD 209458b's transit spectroscopy, 
the vertical extent of the cloud is of great importance. 
It seems likely that any clouds seen during transit spec- 
troscopy will be in a radiative (quiescent) region, in which 
case the cloud particles will likely remain small, in contrast 
to the large sizes (10-200 /im) predicted for convective re- 
gions. Given current uncertainties in cloud modeling (see 
Cooper et al. 2003), much work still needs to be done in 
this area. While Seager & Sasselov (2000) were the first 



to note the importance of clouds in determining the depth 
of transit absorption features, here we note that under- 
standing the vertical distribution of clouds is yet a further 
challenge in EGP transit spectroscopy. 

In this narrow spectral region (581-638 nm), the cross- 
sections have virtually no wavelength dependence. Thus, 
clouds lead to a wavelength-independent absorption of 
stellar light greater than that of the gaseous species alone, 
leading to a uniform increase in the planetary radius across 
much of the wavelength band. In the core of the Na-D fea- 
ture, where the slant optical depth reaches unity at much 
lower pressures, the clouds have no effect. This leads to a 
decrease of the contrast in transit depth between the "in 
feature" and "out of feature" bands. Fig. 3 is analogous 
to Fig. 2, but in this case shows the planet's radius and 
pressure at the terminator as a function of wavelength, 
for our atmosphere including clouds, with the base at 5.4 
mbar. In both cases, the planet's radius at 1 bar has been 
adjusted so that the transit depth remains 1.64%. 

Fig. 4 shows the transit depth versus wavelength for a 
neutral cloud-free atmosphere and a neutral atmosphere 
with the predicted clouds. The cloud-free model gives a 
feature depth of 5.20 x 10^**, deeper than that obtained 
by Hubbard et al. (2001) because of our warmer T-P pro- 
file. The model with clouds decreases this depth signif- 
icantly. With the forsterite cloud base at 5.4 mbar, the 
depth is 2.42 x 10~^, while if the cloud base is at 1.8 
mbar the depth is 0.93 x 10^^. The depth is actually 
smaller than the Charbonneau et al. (2002) observations, 
(2.32 ± 0.57) X 10"^, indicating that if clouds form too 
high in the planet's atmosphere they can dominate over 
the opacity of gaseous species. 

4.2. Sensitivity of Results 

Fig. 1, which shows the T-P profile, along with the pres- 
sure axes on Figs. 2 and 3, highlight how sensitively the 
Na-D feature depth depends on the location of clouds. 
Moving the cloud base ~1 gas scale height can cause an 
appreciable difference in the Na-D feature strength. As 
Fig. 1 and the transit depths for the two models with 
clouds show, a change in the pressure of the cloud base 
of this amount would happen with fairly minor deviations 
from our derived T-P profile for the planet's atmosphere. 
If the cloud base occurs at 1.8 mbar, rather than 5.4 mbar, 
the feature depth moves from comfortably within the error 
bars to significantly below. Even higher clouds (warmer T- 
P profiles) would lead to even weaker Na-D features. At 
several hundred K warmer, no condensate clouds would 
form, as is the case in the model of Barman et al. (2002). 
Perhaps the most important issue is that changing our as- 
sumption of how absorbed stellar radiation is redistributed 
around the planet causes non-negligible changes to the 
planet's T-P profile (e.g., Sudarsky et al. 2002). 

We turn next to the role of ionization, and treat the case 
where the opacity profile varies over the planet. We find 
that the depth to which the Na is ionized changes with 
angle from the planet's subsolar point. In Hubbard et al. 
(2001), Na ionization by UV irradiation was ignored. Here 
we find that ionization of Na probably has a secondary but 
non-negligible effect on the observed Na-D feature. 

5. A SIMPLE NON-LTE SODIUM IONIZATION MODEL 
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The proximity to the parent star subjects the transit- 
ing planet to large amounts of ionizing radiation, making 
equilibrium chemistry treatments inadequate. In partic- 
ular, Na ionization alters the profile of neutral Na. This 
is crucial, since the Na-D line and its transit signature 
are features of the neutral state. A rigorous study of the 
ionization of sodium in the atmosphere of HD 209458b 
would require sophisticated non-LTE and non-equilibrium 
photochemical modeling. However, most of the necessary 
atomic and molecular data (in particidar, collisional rates 
and dissociation parameters) are either poorly known or 
completely unknown. Therefore, we attempt to calculate 
the effects of ionization using simple, parametric models. 

Barman et al. (2002) have performed a detailed study 
of non-LTE effects on excitation and ionization of Na I in 
the atmosphere of HD 209458b. However, they find that 
the UV opacity is dominated by atomic metals (Mg, Al, 
Ca, Fc, and Si) . Because of the much cooler temperatures 
we find, these elements are tied up in condensates and are 
thus depleted from the outer layers of the atmosphere of 
HD 209458b. Therefore, they cannot shield sodium from 
intense stellar ionizing radiation. Moreover, Barman et al. 
(2002) had to parameterize the corresponding collisional 
rates, and they found that free electrons are only due to 
ionization of potassium. 

In view of all uncertainties in the atomic and molecular 
parameters, and especially the UV opacities, we feel that it 
is not warranted to perform a detailed non-LTE study. In- 
stead, we devise a simple, non-LTE analytical model that 
allows us to study the effects of various uncertainties in 
relevant quantities (opacities, free electron densities, col- 
lisional rates) on the degree of ionization of sodium in a 
straightforward way. 

We find that much of the outer layers of the atmosphere 
(P ^1/2 mbar), where Na I has the greatest likelihood 
of having non-LTE level populations, arc ionized by the 
stellar UV flux. We also calculate the change in the ion- 
ization depth as a function of angle from the substellar 
point, which to this point has not been investigated in the 
literature. 

5.1. Description of the Model 

We assume both neutral and ionized sodium are repre- 
sented by one-level atoms, with populations (number den- 
sities) Til and n2- The third unknown of the problem is 
the electron density, rig. The three equations that speciiy 
these unknown number densities are the particle conser- 
vation equation, 

ni+n2 = iVNa , (2) 
charge conservation equation, 

"o = "2 + aNjsSa + l3n2 , (3) 
and the statistical equilibrium equation for sodium, 

ni[Ri2+nc Ci2+nH2 C'u e^] = n2[i?2i+"c C'2i+nH2 C21 eH2] , 

(4) 

where A^Na is the total sodium number density. In Eq. (3), 
we assume that the total number of electrons is given by 
the total number of ionized sodium atoms, plus a contribu- 
tion from an unspecified source that is given through some 
fraction a of the total number of sodium atoms (that is, it 
scales as the total density), plus another empirical contri- 
bution that scales as ionized sodium, with a multiplicative 



factor (3. Setting a = /3 = stipulates that free electrons 
are provided only by ionization of sodium. 

In the statistical equilibrium equation (4), R12 and R21 
are the photoionization and radiative recombination rates, 
respectively; ^10^^12 and ncC'21 are the collisional ionization 
and recombination (three-body recombination) rates with 
free electrons. Since the dominant species in the atmo- 
sphere is H2, we have to consider ionization of sodium by 
collisions with H2 and its inverse process as well. In the 
absence of any knowledge of relevant rates, wc assume, 
similar to Barman et al. (2002), that the collisional rate 
is given by a scaled electron rate with some factor eH2- 
It is reasonable to consider this factor to be the ratio of 
reduced masses of the Na -I- H2 pair to that of Na + elec- 
tron, which is about 1/60. However, we keep eH2 as a 
free parameter, allowing us to study the effects of the H2 
collision efficiency on Na ionization. 

By detailed balancing arguments (e.g., Mihalas 1978), 
the inverse rates arc given by 

C21 = Ci2(ni/n2)* = ne0(T)Ci2, (5) 

where (ni/n2)* denotes the LTE value, and (j) is the Saha- 
Boltzmann factor, given by 

<I>{T) = 2.0706 X 10-16 (ffi/52) exp(/ii/o/fcT) , (6) 

where T is the temperature, gi and 92 the statistical 
weights ((?i = 2,(72 = 1 for sodium), the ionization 
frequency, and h and k the Planck and Boltzmann con- 
stants. 

The radiative rates are given by 



-R12 



- ■ 



ail.) 



Ji,dv , 



(7) 



and 



R2i=ne(l>{T)^ / 



( + j exp{-hu/kT)di' , 

(8) 

where a{u) is the photoionization cross-section, and J^, is 
the mean intensity of radiation. However, Eq. (8) takes 

into account only recombinations to the ground state of 
Na I. In reality, recombinations can occur to any state, 
and at low temperatures these processes are very efficient. 
Wc therefore take the recombination rate to be 



C21 = 5.461 X lO-^y\^/T/n{l + Vr/n) 



V 0.825 



(9) 



as given by Verner & Ferland (1996). 

The radiation intensity should in principle be deter- 
mined by solving the radiative transfer equation. We could 
do so, but in view of gross uncertainties in relevant UV 
opacities, this is not warranted. Instead, we simply as- 
sume that the specific intensity is given as a sum of at- 
tenuated stellar irradiation intensity plus a thermal com- 
ponent given by the Planck function at the local temper- 
ature. (This part is negligible at surface layers of interest, 
but gives the correct behavior of intensity at depth.) To al- 
low for uncertainties in UV opacities, we introduce another 
multiplicative factor, Wadd) which represents an additional 
attenuation of radiation, perhaps by an unspecified opac- 
ity source. 

The intensity as a function of depth is given by 



Lf,{.T^) = Wadd 4^ exp(-T^/i) + B^{T) , 



(10) 
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where ^ is the cosine of the angle with respect to the sur- 
face normal, t^(P) is the optical depth (in the normal 
direction), and PJ^^^ is the specific intensity of the stellar 
radiation at the top of the planetary atmosphere. We use 
a Kurucz model for a GOV star with a dilution factor of 
m; = 1.4 X 10~^ due to the planet-star distance. For sim- 
plicity, we assume isotropic irradiation, which gives 

J- = lj^ I^^^^d^i = ^ w^d IZM-r) + B^T) , (11) 

where E2 is the second exponential integral. 

The electron collisional ionization rate is assumed to be 
given by Seaton's formula, 

C12 = 1.55 X 10^^T-^^^ga{iyo) exp(-Mo)/«o , (12) 

where (t(i'o) is the photoionization cross-section at the 
edge, g = 0.1, and uq = hvo/kT. 

5.2. UV Opacities 

Since we find no free heavy elements in atomic form 
above the clouds in which these elements are sequestered, 
any opacity that may affect Na ionization must be pro- 
vided by molecular species. We have searched the liter- 
ature, and found only two equilibrium species that con- 
tribute significantly to the opacity in the near UV H2O 
(Yoshino et al. 1996), and H2S (Lee et al. 1987). Raylcigh 
scattering on H2 is much less important at this range. For 
both H2O and H2S data are only avaivable at 295 K. We 
plot in Fig. 5 the total opacity, plus the individual contri- 
butions of H2O, H2S, H2, and a maximum possible con- 
tribution of Na I (that is, assmning that sodium is com- 
pletely neutral) at upper layers of the atmosphere of HD 
209458b, as they follow from the atmosphere model. The 
corresponding Na I photoionization cross-section is taken 
from the Opacity Project database (Cunto et al. 1993). 

5.3. Representative Results 

The basic set of equations (2) - (4), together with aux- 
iliary expressions (5) - (11), can be cast into a single cu- 
bic equation for n2/N^SL- All the necessary quantities, T, 
A^Na, '^■(112), n(H20), and n{H2S), as functions of depth, 
are provided by our model atmosphere for HD 209458b. 

Our nominal model, without any empirical modifica- 
tions, uses a. = f3 = 0, £^12 = 1/60, and Wadd = 1- The 
depth dependence of Na ionization for this model is shown 
in Fig. 6. This calculation corresponds to the planet's sub- 
solar point. For the nominal model, the pressure at which 
the number densities of Na I and Na+ are equal is 1.6 
mbar. At greater pressures, neutral Na dominates. Deep 
in the planet's atmosphere, at pressures > 100 bar, ion- 
ization of Na again becomes appreciable due to the higher 
temperatures. Also shown is the ionization depth depen- 
dence without shielding due to the other UV absorbers. 
The pressure at which Na is 50% ionized increases to 3 
mbar. 

We next test our predictions for Na ionization depth by 

varying several of the parameters mentioned earlier. These 
trials can be seen in Fig. 7. In the first case a — I, which 
represents the case where there is another donor, with a 
total abundance equal to that of sodium, which is fully 
ionized. The model with /3 = 1 represents an additional 



empirical donor of free electrons that has the same abun- 
dance as ionized sodium. The model with Wadd = 0.1 rep- 
resents the situation where unknown UV absorbers leave 
only 10% of the calculated incoming ionizing flux available 
to ionize Na. The electron donor factors seem like reason- 
able upper limits because species more abundant than Na 
in the planet's atmosphere tend to have 2-3 times larger 
ionization energies and consequently are difficult to ion- 
ize. These factors have a modest effect on the pressure 
at which Na is 50% ionized. Changing the value of eH2 
(which characterizes the H2 collision efficiency) by factors 
of even 100 had negligible effects on the ionization depth 
for pressures < 100 mbar. 

5.4. Ionization on the Planet's Limb 

The important drawback of the ionization model out- 
lined above is that the calculations are representative of 
the subsolar point only. We are interested in ionization 
at the planet's limb. At a given location in the planet's 
atmosphere an angle from the subsolar point, the ab- 
sorbed radiation of a given surface area will be diluted by 
a factor cos0, meaning the ionization depth will not be as 
great. Even this geometric view is too simple, as it would 
imply no photoionization at the planet's terminator. To 
overcome this shortcoming, we have employed a somewhat 
different ionization model that can predict the ionization 
depth vs. angle around the planet. 

In this angle-dependent 2-D ionization model (which we 
will call model 2, while the prior model is model 1), wc as- 
sume incoming photons with energy greater than 5.14 eV 
(the ionization energy of Na) and less than 11.2 eV, the 
energy at which significant H2 absorption begins, ionize 
atomic Na. Na"*" is the only ionized species and conse- 
quently the number of free electrons is equal to the num- 
ber of Na"*" ions. We use 2- and 3-body radiative recom- 
bination coefficients from Verner Ferland (1996) and 
Gudzenko et al. (1974). Incoming photons at some oblique 
angle relative to the surface can ionize all the Na along a 
chord until, at some pressure, most ionizing photons (per 
second) have been used up, and recombinations (per sec- 
ond) can match ionizations (per second). This is the ion- 
ization depth. Much as for a Stromgren sphere, we find 
a sharp boundary between regions of the planet's atmo- 
sphere with complete ionization and with no ionization. 
(Calculations from model 1 find a more gradual onset of 
neutral Na with depth.) We calculate the depth to which 
the atmosphere is ionized at several hundred angular po- 
sitions from the planet's substellar point. We find that Na 
is ionized to a depth of ~12 mbar at the planet's substellar 
point (this was 1.6 mbar for model 1 with shielding due 
to UV absorbers, and 3 mbar without), and as one moves 
away from this point, the pressure to which Na is ionized 
decreases. At the terminator (90 degrees from the subsolar 
point) the ionization depth is 1.7 mbar. Fig. 8 shows rep- 
resentative calculations for the change in ionization depth 
with angle near the terminator. At an angle of ~10 degrees 
past the terminator, no Na is ionized. The model possesses 
cylindrical symmetry about the planet's subsolar- antisolar 
axis. The removal of the opacity of atomic Na in the 10 
mbar to 1 /ubar region of the atmosphere produces a neg- 
ligible difference in the T-P profile. 

The column density of atmosphere through which pho- 
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tons travel until one reaches the ionization depth of Na 
can be calculated at the subsolar point for model 1, and 
at any angle, including the subsolar point for model 2. Us- 
ing the ratio of the column depths of the two models at 
the subsolar point, one can scale the results from model 
2 to calculate the ionization depth around the planet for 
model 1. Fig. 8 shows the ionization depth versus angle 
for the scaled model 1 results. 

Fig. 9 depicts how ionization leads to a shallower tran- 
sit depth. The differences between the Charbonneau et al. 
(2002) "in feature" and "out of feature" bands are, for 
model 1, 4.72 x 10"^ and for model 2, 4.83 x lO"*. Nei- 
ther of these models include the predicted clouds, so that 
the relatively minor effect due to ionization can be shown 
clearly. Fig. 8 shows a model in which we increased the 
pressure level to which Na is ionized, in order to gauge 
how much ionization is needed to match the Charbonneau 
et al. (2002) observations without the opacity of clouds. 
The feature depth for this model is 2.87 x 10~^. 

6. CONCLUSIONS 

We find that HD 209458b's atmosphere is cool enough 
for most heavy elements to have been incorporated into 
condensates but still warm enough for Na to remain in 
atomic form. As seen in Fig. 1, there is a range of ~600 K 
in the mbar region of the planet's atmosphere where this 
situation holds. High clouds arise naturally in the planet's 
atmosphere near ^1-5 mbar. Although the abundances of 
condensable vapor of the chief cloud-forming species are 
small, the particles that form in this radiative region are 
small and have large scattering cross-sections. Clouds lead 
to an increase in the apparent planet radius outside of the 
Na-D core, and on their own they could easily explain the 
observations. In addition, clouds can be so high that the 
predicted sodium feature is less than observed. 

With our sinipk^ non-LTE model, we have shown that 
stellar UV flux will ionize Na to ~l/2 mbar on the limb 
in the planet's atmosphere, reducing the strength of the 
589-nm absorption feature. The change in the Na-D fea- 
ture due to ionization is fairly insensitive to the derived 
T-P profile. Ionization is a secondary effect, and it alone 
cannot explain the weak Na-D features. This conclusion 
was also reached by Charbonneau et al. (2002) in their 



analysis. Since the majority of Na on the planet's limb 
at pressures ^1/2 mbar is ionized by stellar UV flux, the 
possible effects of neutral Na non-LTE level populations 
found by Barman et al. (2002) are reduced. In addition, 
since the T-P profile we derive is cooler, a region where 
Na LTE level populations may not hold will be pushed out 
to lower pressure. 

For our T-P profile, we find an opaque forsterite cloud 
base at a pressure of several mbar. Predicting the exact 
location of the cloud base and vertical distribution of cloud 
material is challenging, but the Na-D feature is quite sen- 
stive to these factors. An atmosphere incorporating the 
forsterite cloud at 1.8 mbar and ionization gives a Na-D 
feature less than observed (due to clouds, not ionization), 
but for a cloud base of at 5.4 mbar, we obtain a feature 
depth that matches observations with or without ioniza- 
tion. 

More advanced modeling of photochemistry may reveal 
other non-equilibrium compounds, perhaps including a 
photochemical haze, which could both increase the impor- 
tance of particulate absorption and decrease the extent to 
which the atmosphere is ionized. The potassium feature at 
770 nm should be more sensitive to ionization than is Na, 
because the former's ionization energy is 0.8 eV less than 
that of Na. We are currently updating our UV opacity 
database and collecting all available collsional parameters. 
Once this is completed, it will be worthwhile to revisit the 
sodium ionization problem, taking into account many lev- 
els of Na I, and computing consistent non-LTE level pop- 
ulations of sodium together with the radiation field. More 
advanced global circulation models for HD 209458b will 
be of great interest, as they will give a better indication of 
the temperatures expected on the limb of the planet. It is 
likely that observations will push modeling in new direc- 
tions, and data obtained across many wavelengths, both 
from the ground and from space, including predicted K, 
H2O, and CO features, will help gauge the importance of 
processes at work in HD 209458b's atmosphere and decide 
between explanations of the Na-D feature depth. 

The authors acknowledge the valuable insight gained 
through discussions with Don Hunten and John Plane and 
support from grants NAG5-10760, -10450, and -10629. 



REFERENCES 



Barman, T. S., Hauschildt, P. H., Schweitzer, A., Stancil, P. C, 

Baron, E., & AUard, F. 2002, ApJ, 569, L51 
Brown, T. M. 2001, ApJ, 553, 1006 

Brown, T. M., Charbonneau, D., GiUiland, R. L., Noyes, R. W., & 

Burrows, A. 2001, ApJ, 552, 699 
Burrows, A., Hubbard, W. B., Lunine, J. I., & Liebert, J. 2001, 

Reviews of Modern Physics, 73, 719 
Castor, J. I., Dykema, P. G., & Klein, R. I. 1992, ApJ, 387, 561 
Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, M. 2000, 

ApJ, 529, L45 

Charbonneau, D., Brown, T. M., Noyes, R. W., & GiUiland, R. L. 

2002, ApJ, 568, 377 
Cooper, C. S., Sudarsky, D., Milsom, J. A., Lunine, J. I., & Burrows, 

A. 2003, ApJ in press, astro-ph/0205192 
Cunto, W., Mendoza, C., Ochsenbein, F., & Zeippen, C. J. 1993, 

A&A, 275, L5+ 

Gudzenko, L. I., Evstigneev, V. V., Filippov, S. S., & Yakovlenko, 

S. L 1974, Teplofizika Vysokikh Temperatur, 964 
Guillot, T. & Showman, A. P. 2002, A&A, 385, 156 
Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, ApJ, 

529, L41 

Hubbard, W. B., Fortney, J. J., Lunine, J. I., Burrows, A., Sudarsky, 
D., & Pinto, P. 2001, ApJ, 560, 413 



Hubeny, I. & Lanz, T. 1995, ApJ, 439, 875 

Lee, L. C, Wang, X., & Suto, M. 1987, J. Chem. Phys., 86, 4353 
Mazeh, T., Naef, D., Torres, G., Latham, D. W., Mayor, M., Beuzit, 

J., Brown, T. M., Buchliavc, L., Burnet, M., Carney, B. W., 

Charbonneau, D., Drukier, G. A., Laird, J. B., Pepe, F., Perrier, 

C, Queloz, D., Santos, N. C, Sivan, J., Udry, S., & Zucker, S. 

2000, ApJ, 532, L55 
Mihalas, D. 1978, Stellar atmospheres 2nd edition (San Francisco, 

W. H. Freeman and Co., 1978. 650 p.) 
Richardson, L. J., Deming, D., Wiedemann, G., Goukenleuque, C, 

Steyert, D., Harrington, J., & Esposito, L. W. 2002, in ApJ in 

press, astro-ph/0210612 
Seager, S. & Sasselov, D. D. 2000, ApJ, 537, 916 
Showman, A. P. & Guillot, T. 2002, A&A, 385, 166 
Sudarsky, D., Burrows, A., & Hubeny, L 2002, ApJ submitted, astro- 

ph/0210216 

Sudarsky, D., Burrows, A., Pinto, P. 2000. ApJ, 538, 885 

Verner, D. A. & Ferland, G. ,1. 1996, ApJS, 103, 467 

Yoshino, K., Esmond, J. R., Parkinson, W. H., Ito, K., & Matsui, T. 

1996, Chemical Physics 211 (1996): 387; also in Chemical Physics 

215 (1997): 429., 211, 387 



On the Detection of Sodium in HD 209458b 



7 




500 1000 1500 2000 2500 

T (K) 

Fig. 1. — Temperature-pressure profile for the atmosphere of HD 209458b (sohd line). The dotted line is the T-P profile we derive without 
the opacity due to clouds. The condensation curves for Na2S, forsterite {Mg2Si04), and iron are shown as dashed lines. When cloud opacity 
is taken into account the forsterite cloud base moves from 5.4 mbar to 1.8 mbar. There is a wide range of possible T-P profiles for which 
forsterite and iron clouds form high in the planet's atmosphere but Na remains in atomic form. We used a surface gravity of 980 cm/s^ for 
the planet. 



8 



Fortney, et al. 



1 .400 
1 .390 

1 .380 

1 .370 



^ 1.360 
o 



1 .350 ^. 

1 .340 

1 .330 
580 




600 620 

A fnm 



640 



Fig. 2. — Planetary radius versus wavelength, also showing the pressure at the planet's limb corresponding to the derived radius, for a 
neutral cloudless atmosphere. (1 Rj = 71492 km, Jupiter's equatorial radius at 1 bar.) The radius values are the radii derived from calculating 
the apparent cross-sectional area of the spherical planet, and this radius is essentially identical to the radius at a slant opactical depth of 2/3. 
The wavelengths shown encompass the Brown et al. (2001) observations. During the transit, an observer is sensitive to pressures from tens 
of mbar to ~10 /ibar. The average pressure across this wavelength band is 15 mbar. 
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Fig. 3. — Similar to Fig. 2, but for a neutral atmosphere with a forsterite cloud base of 5.4 mbar. The pressure axis changes due to the 
different opacity in the planet's atmosphere, as compared to Fig. 2. The clouds lead to a wavelength-independent increase in planetary radius. 
The average pressure across this wavelength band is 1 mbar. 
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Fig. 4. — Transit depth (at midtransit) versus wavelength for several models. Curves appear oflFset because the transit depth is required 
to be a fractional change in flux of -0.0164 from 581 to 638 nm, the full Brown et al. (2001) observation band. The two black plusses mark 
the wavelength extent of the "in feature" band from Charbonneau et al. (2002), while the wavelengths between the plusses and diamonds, 
on either side of the "in feature" band, make up the "out of feature" band. The observed difference in transit depth between these bands is 
(2.32 ± 0.57) X 10~^. Shown are neutral models with and without clouds. The Na feature depth for neutral cloudless model (thin solid line) 
is 5.20 X 10~*. For a forsterite cloud base at 5.4 mbar, the feature depth is 2.42 x 10~* (dashed line), and 0.93 x 10"'' for a cloud base at 
1.8 mbar (thick black line). 
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Fig. 5.— 



Opacity at the outer layers of the atmosphere of HD 209458b. 
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Fig. 6. — Influence of shielding on the Na ionization depth due to UV absorption by H2O, H2S, and Rayleigh scattering. The thick solid 
line is our nominal model, while the dashed line is the same as the nominal model, but with Na able to absorb all ionizing radiation. The 
thin horizontal dotted line shows the location of 50% Na ionization. 
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Fig. 7. — Influence on the ionization depth of factors a, (3, Waddi leaving other parameters at their default values. As described in the text, 
these factors are meant to show the sensitivity of our results to possible unknown additional free electrons or additional UV opacity. The 
thick solid curve is our default calculated ionization depth curve, {a = f3 = eH2 = 0, tUadd = 1) The thin solid curves are for a = 1 and /? = 1, 
while the dashed curve is TOadd = 0.1. The thin horizontal dotted line shows the location of 50% Na ionization. All these factors would push 
the ionization depth out to lower pressures. 
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Fig. 8. — Sodium ionization depth versus angle from the planet's substellar point. The dotted vertical line at 90 degrees indicates the 
terminator. The dashed curve is for model 1, while the solid curve is for model 2. All Na is ionized to the left of this curve and neutral to 
the right of the curve. The dot-dash curve shows an increased ionization depth at every angle, an example of an ionization depth curve that 
would be necessary to match the observations in the absence of any cloud opacity. 
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Fig. 9. — Transit depth (at midtransit) versus wavelength for several models, similar to Fig. 4. The Na feature depth for neutral eloudless 
model (thin solid line) is 5.20 X 10^^. For the ionized cases, the feature depth is 4.72 X 10^'* for model 1 (dashed line), and 4.83 X 10~* for 
model 2 (thiek solid line), respeetively. Even though the ionization depth is greater for model 2 (see Fig. 8 and the text), the "in" - "out" 
feature difference is greater, because in model 2, the Na wings are affected, but in model 1 they are not. 



